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Abstract 

We construct a gauge-mediated SUSY breaking model with vector-like matters 
combined with the Peccei-Quinn mechanism to solve the strong CP problem. The 
Peccei-Quinn symmetry plays an essential role for generating sizable masses for the 
vector-like matters and the //-term without introducing dangerous CP angle. The 
model naturally explains both the 125GeV Higgs mass and the muon anomalous 
magnetic moment. The stabilization of the Peccei-Quinn scalar and the cosmology 
of the saxion and axino are also discussed. 



1 Introduction 



Recently ATLAS and CMS collaborations have discovered the Higgs boson at the mass 
around 125 GeV [U[2]. In conventional gauge mediated supersymmetry (SUSY) breaking 
(GMSB) models |3j, it is difficult to explain the Higgs mass of 124-126 GeV, unless the 
SUSY particles are as heavy as 10-100 TeV. Obviously such heavy SUSY particles are not 
favored from a viewpoint of naturalness. Moreover, there is another possible indication of 
TeV scale SUSY particles : the muon anomalous magnetic moment (muon g — 2). In fact 
the experimental value of the muon (7 — 2 is deviated from the SM prediction at about 3 
u level [4j [5] . This deviation can be naturally explained with TeV scale SUSY particles 
and relatively large tan f3. 

One of the easiest way to raise the Higgs mass up to 125 GeV while keeping the SUSY 
particle masses around 1 TeV is to introduce additional vector-like matter [6j , which couple 
to the minimal SUSY standard model (MSSM) Higgs via yukawa interactions. Similarly 
to the top yukawa coupling which radiatively increases the Higgs mass, the new yukawa 
coupling gives an additional potential for the Higgs fields and the Higgs can be heavier 
than the MSSM case. In order to avoid the gauge anomaly, the vector-like matter may be 
in the 10 and 10 representations of SU(5) grand unified theory (GUT) gauge group. This 
type of extension of the MSSM has been recently discussed in Refs. [TJ El El ED] in the light 
of recent ATLAS and CMS results. In particular, it was pointed out that the 125 GeV 
Higgs and the muon g — 2 can be explained simultaneously in this class of models [SI EI] • 

From a viewpoint of model building, however, this is far from complete. First, the 
/i/-B/i-problem in GMSB must be solved in order to obtain a correct electroweak sym- 
metry breaking (EWSB) minimum. We cannot discuss the SUSY CP problem unless the 
mechanism for generating fi/Bfi is specified!*] Second, the SUSY masses of vector-like 
matter, which can in principle take arbitrary values, must also happen to be around the 
weak scale in order to raise the Higgs mass. Moreover, in general there can be both 
the up-type and down-type Higgs couplings to the vector-like matter, the latter of which 
tends to decrease the Higgs mass. One must somehow tune the latter coupling so as not 



1 In the gravity-mediation models, /1 and B[i terms as well as the SUSY masses of the vector-like 
matter can be generated [7] by the Giudice-Masiero mechanism [11] . 
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to affect the Higgs potential. Finally, the strong CP problem was not addressed in these 
frameworks. 

In this paper we consider the extended GMSB model with Peccei-Quinn (PQ) symme- 
try [12] in order to deal with above mentioned problems. MSSM fields as well as vector- 
like matter are charged under the U(1)pq, so that the /x-term and the mass terms for 
the vector-like matter are forbidden. Then they are generated by the vacuum expectation 
value (VEV) of the PQ scalar field at a correct scale for the PQ scale of 10 9 -10 n GeV |13j . 
The spontaneous PQ symmetry breaking predicts an almost massless Nambu-Goldstone 
boson, called axion, which dynamically solves the strong CP problem [T2l [Ti] . Dangerous 
CP angle is not introduced, hence it also solves the /x/-B/x-problem and the SUSY CP 
problem. Unwanted couplings between down-type Higgs and the vector-like matter are 
forbidden by the PQ symmetry. 

This paper is organized as follows. In Sec. [2] our model is introduced, and it is shown 
that the Higgs mass of 125 GeV and the muon g — 2 can be explained simultaneously. In 
Sec. [3] a mechanism for stabilizing the PQ scalar is described. Cosmological constraints 
on our model are also discussed. Sec. H] is devoted for summary and conclusions. 

2 The extended GMSB model with PQ symmetry 

We consider a model of gauge-mediated SUSY breaking with global U(1)pq symmetry. 
The superpotential of the model is given by 

W = Wmssm + + W PQ+ 

mess; V / 

where Wmssm contains Yukawa interaction terms in the MSSM (except for the /x-term), 
Wpq +mcss is the superpotential for the PQ and messenger sector, and W ext consists of the 
additional vector-like matter. They are given by 

W ext = ^f-H^Hi + -^-(A 2 Q'Q' + X 3 U'U' + X 4 E'E') + Y'Q'H U U', (2) 

IVlp IVlp 

and 

WpQ+mess = &0# P q^pq + /tX* mess ^ mess , (3) 
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Table 1: Charge assignments on chiral superfields fields in the model under the £/(1)pq 
and R-parity (+ : even, — : odd). 5m and 10m ar e the MSSM matter fields. 



where is a PQ symmetry breaking field and (Q', U', E') and (Q', U', E') are the extra 
vector-like matter, which transform 10 and 10 under the £77(5) GUT gauge group, re- 
spectively, and M P is the reduced Planck scale. The PQ quarks, \&pq and ^pq, and the 
messenger fields, ^ mess and ^mcss transform as 5 and 5 under the SU(5) GUT. The spurion 
field X gives the SUSY breaking mass to the messenger. The PQ scalar obtains a VEV 
of order of 10 10 — 10 12 GeV as explained in the next section. Thus it spontaneously breaks 
the U(1)pq symmetry and the associated NG boson behaves as the axion which solves the 
strong CP problem [TBI [16]. Then the SUSY masses for the Higgs and vector-like matter 
are generated as fi ~ M<yp<,E? ~ A 4 (0) 2 /Mp for (0) ~ 10 10 - 10 12 GeV. The PQ charge 
assignments on these fields are summarized in Table. [U 

Couplings between the Higgs and the extra matters give large radiative corrections to 
the lightest Higgs mass. The last term in Eq. (j5J), Y'Q'H U U' , gives positive corrections to 
the Higgs mass squared as [TTJ, [18] § 

Am\ - 4 sin 4 bL M _ I U - MV 5 - ^] 

where Mp(~ Mqijji) is a fermionic mass of the extra matter and Ms is the average of the 
scalar masses, M| = Mq, w + m 2 oft where m soft denotes the SUSY braking contribution. 
Here we consider the decoupling limit. On the other hand, if the coupling like Y"Q'HdU' 
would exist, it gives negative contributions to the Higgs mass squared, which can be 

2 The finite corrections, which arise through trilinear couplings, A'Y' Q' H U U' are not shown here. In 
numerical calculations, they are included. 



(4) 
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significantly large for Y" ~ 1 and large [i parameter. The correction is given by [HI [E] 

where corrections suppressed by tan/3 are neglected. Moreover, a CP violating phase 
exists, since phases of Mqi, Mjjr, Y' and Y" can not be removed simultaneously. In our 
model, Y" is of the order of (0) 6 /Mp thanks to the PQ symmetry and hence highly 
suppressed. Therefore this effect is negligibly small. Note that there are no significant 
corrections to couplings, h-gg and fo-77, as long as Y" is negligibly small; the corrections 
are approximately proportional to A = J2 i= i 2 (d log Mp ^/(d log v u ) (see e.g., [19]), where 
Mp,i is the eigenvalue of the fermion mass matrix of the extra-matters, and A vanishes 
for Y" = 0. 

Due to the additional contributions, the Higgs mass (m^) can be easily around 125 GeV 
with neither the enhanced trilinear coupling of the stop nor heavy stops. This is welcome 
in terms of the muon g — 2. In fact, the experimental value of the muon g — 2 is deviated 
from SM prediction at 3.2<r [I]: 

fl ex P _ a SM = (26.1 ± 8.0) x 1(T 10 . (6) 

The deviation can be naturally explained for large tan and relatively small soft SUSY 
breaking mass. In the MSSM without large A t , rrih — 125 GeV requires the stop mass 
of (9(10) TeV. Consequently, the Higgs mass and the muon g — 2 can not be explained 
simultaneously. Remarkably, in the model with the vector-like matter, the deviation 
can be explained consistently with the Higgs mass m>h — 125 GeV in GMSB due to the 
additional contributions to the Higgs mass jSJUU]. In the numerical calculation, the SUSY 
mass spectrum is calculated by Suspect package [20], which is modified to include 2- loop 
renormalization group equations for the extra matter. The Higgs mass and the muon 
g — 2 are evaluated by FeynHiggs package [21] . 

In Fig. [TJ contours of the Higgs mass and muon g — 2 for different messenger scales 
are shown. The gray regions are excluded by the constraint from the charge breaking 
vacuum [22] (see Ref. [23] for GMSB with vector-like matter) or/and the LEP bound on 
the stau mass, rrif < 87.4 GeV [24J. Due to additional negative corrections to m 2 H from 
the vector-like matter, the predicted value of \i parameter tends to be larger than that of 
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MSSM. As a result, a trilinear coupling of the stau, 



2M W 



H tan /Stlt^/i + h.c, 



(7) 



becomes large and the charge breaking minimum, which can be deeper than the elec- 
troweak symmetry breaking minimum, might be generated. 

The Higgs mass is predicted to be 124 — 126 GeV in the red bands. We have taken 
the SUSY masses for the vector-like matter as Mq> = Mjj' = 600 (1200) GeV in the left 
(right) bands. In the orange (yellow) region, muon g — 2 is explained at 1 a (2 a) level 
and the corresponding gluino mass is m g < 1.3 (1.8) TeV. Since the current results from 
the SUSY searches exclude the region with rrig < 1 TeV as discussed later, the messenger 
scale should be lower than ~ 10 7 GeV (~ 10 10 GeV for 2 a) if the SUSY is responsible 
for the muon g — 2 anomaly. On the blue dashed line, the lightest nuetralino and the 
stau are degenerate in the mass. The region above (below) the line, NLSP is the stau 
(neutralino). 

The line B(M mess ) = corresponds to the vanishing 5-term at the messenger scale; 
the Higgs S-term is generated radiatively through the gaugino masses, and the successful 
electroweak symmetry breaking is achieved for relatively large tan /3 [25j. In this case the 
phases of the B-term and gaugino mass are aligned, and there is no SUSY CP problem 
as well as fi/Bfi problem. 

The contours of the Higgs I?-term at the messenger scale, B(M mess ), are shown in 
Fig. [2j As the messenger scale becomes high, the radiative correction from the renor- 
malization group evolution between the messenger scale and the electroweak symmetry 
breaking scale also becomes large. The successful electroweak symmetry breaking requires 
the following relation at the weak scale 



where the radiative corrections are neglected. For large tan/3, the right hand side of the 
above equation should be small; the Higgs .B-term should be small for fixed values of fi 2 
and m 2 rr . As a result, the low messenger scale is favored, otherwise the generated 5-term 



is too large to satisfy the stationary condition for the successful electroweak symmetry 
breaking. 



tan/3 Bjj, 



(8) 



l + tan 2 /3 2 1 /x | 2 + m? H + m 
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M mcss = 4 x 10 s GeV M mess = 10 6 GeV 
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Figure 1: The contours of the Higgs mass and muon g — 2 on m^-tan (3 plane for different 
messenger scales. The horizontal axis displays the gluino mass in the unit of GeV. On 
the blue dashed line, the lightest neutralino mass and the lightest stau mass are equal, 
mfj = Tn Xo . In the red bands, the Higgs mass of 124 GeV < rrih < 126 GeV is explained 
for Mq^jji = 600 GeV and 1200 GeV. The gray region is excluded by the vacuum stability 
bound and/or the bound from OPAL experiment [24|. 
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However even in the case of B(M mcss ) ^ 0, the successful electroweak symmetry 
breaking is realized without generating a dangerous CP violating phase. The additional 
contribution to the Higgs B-term arises through the following interaction, 

g D Q'H d y D + g L E'H d V L , (9) 

or 

g D Q'H u y D + g L E'H d y L , (10) 

where ^ ^ and are the parts of the SU(5) multiplet \l/ mcss , transforming 3* x 1 and 
1x2 under SU (3)c x SU (2) l, respectively. The above interaction (J9]) ( (JTOjl ) is allowed by 
choosing the PQ charge of the messenger as g meS s = — 2 — 3g#/2 (g mess = 2 — 3g#/2). Since 
the phases of the couplings gjj, gi and k (in Eq.(J3])) can be removed simultaneously, the 
above interaction does not generate any CP violating phases. The additional contribution 
to the Higgs .B-term is given by c 

6B(M mess ) ~ -J^(3gD + ^!)Amess, (11) 

with A mcss = Fx/X, where we have shown the only leading contribution. Since SB is 
negative, the region above the line B(M mcss ) = in Fig. [2] is also consistent with the 
correct electroweak symmetry breaking. 

Finally, let us comment on the constraints from the recent results of the SUSY searches. 
In our model, the muon g — 2 and the Higgs mass can be explained simultaneously in wide 
range of the messenger scale, M mess ~ (a few) x 10 5 GeV - 10 10 GeV. For the messenger 
scale as low as M mess ~ 10 5 GeV, the gravitino mass can be as light as m 3 / 2 ~ 10 eV if 
the dominant SUSY breaking comes from the F-term of X {Fx). In this case, the next 
to the lightest SUSY particle (NLSP) decays into gravitino and SM particle promptly 
(ct nlsp < 0.1 mm). When the stau is the NLSP, the typical SUSY signal contains 
large missing transverse momentum, jets and r leptons. In such a case, the current 
bound of the gluino mass is obtained as rrig > 1 TeV |2B]. On the other hand, In the 
region where the bino-like NLSP decaying into the photon and the gravitino, the signal 
contains large missing transverse momentum and photons, and the constraint is obtained 



Soft mass squared of Hd(H u ) is also modified by the interaction ©(dUO- 
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Figure 2: The contours of the Higgs 5-term at the messenger scale, -B(mess), is shown in 
the unit of GeV on M mess -tan (3 plane. 



as nig > 1.2 TeV for nig > 1.6 TeV [27J. All regions consistent with muon g — 2 at the la 
level are already excluded. 

On the other hand, if the messenger scale is high enough and/or the SUSY breaking 
F-term (F tota i) is much larger than Fx, the decay length of the NLSP is longer than 
the detector size. Therefore the NLSP can be regarded as the stable particle inside the 
detector. The stau NLSP region which is consistent with muon (7 — 2 is expected to be 
excluded [TU]. When the neutralino is the NLSP, the typical signal contains multi-jets 
and missing transverse momentum as in the case of mSUGRA models. In such a case, 
the current bound is estimated as m g > 1 TeV for rrig > 1.4 TeV [28J (see also Ref. [10J). 

In Fig. [31 contours of the squark masses are shown on m g — M mess plane. Current 
bounds from LHC SUSY search can be avoided for rrig > ITeV and m g > 1.5 TeV. On 
the other hand, the muon g — 2 at the la (2a) level and the 125 GeV Higgs boson are 
explained for rrig < 1.3 (1.8) TeV. In our model, the gluino mass m g > 1 TeV corresponds 
to rrig > 1.6 TeV depending on the messenger scale. Therefore, the Higgs mass and muon 
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g — 2 are explained simultaneously, while satisfying the current bound from the LHC 
SUSY search. 

In summary, our model has following properties. 

• The size of /i-term and the masses of vector-like matter are controlled by the PQ 
symmetry. They are naturally of the order of TeV for the phenomenologically viable 
PQ scale of 10 9 -10 12 GeV, hence solves the /z-problem. The B-term can be either 
zero at the messenger scale or may be generated through the interaction ([9]) or ( fill . 

• Due to the radiative correction from additional vector-like matter with masses of 
(D(TeV), the lightest Higgs mass can easily be as heavy as 125 GeV. The unwanted 
coupling of the vector-like matter to the down type Higgs is forbidden by the PQ 
symmetry. The muon anomalous magnetic moment can also be explained simulta- 
neously while current SUSY searches at the LHC can be avoided. 

Here are additional comments. 

• The PQ symmetry is anomaly- free for a particle content included in Wmssm and W ext . 
Thus we have introduced one pair of PQ quarks, ^pq and ^pq, which transform as 
5 and 5 under SU(5) in order to make the U(1)pq anomalous. This solves the strong 
CP problem. The domain wall number is equal to one and hence we do not suffer 
from the cosmological domain wall problem. It also stabilizes the PQ scalar at the 
scale of 10 9 -10 12 GeV, as will be seen in the next section. 

• The perturbativity of the gauge couplings is maintained if the mass of heavy quarks, 
Mp Q = k{<f>), is larger than about 10 10 GeV for iV mcss = 1 and M mess = 10 6 GeV. 
Larger M mess leads to the looser constraint on Mpq. 

• The seesaw mechanism can work for = —4/5. The extra matter can decay into 
SM particles through higher dimensional operators, with the life time of ~ 1 sec, and 
it can be consistent with the successful prediction of the big-bang nucleosynthesis 
(BBN) (see Appendix B). 

• The extra-matters can have odd R-parity. In this case, the extra-matters mix with 
SM particles, and thus they decay into SM particles promptly, (see also recent 
discussions about discovery of the extra- matters [TOj, 29J) 
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Figure 3: The contours of the squark mass on mg-M mess plane. The squark masses are 
shown in the unit of GeV. 
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3 Stabilization of the PQ scalar and cosmology 
3.1 Stabilizing the PQ scalar 

In this section we discuss how to stabilize the PQ scaler at an appropriate scale. The 
stabilization of the PQ scalar in the framework of GMSB was discussed e.g. in Refs. J30J 
EH [321 [331 EH [35] . For k\<f>\ < M mess (= kX), the Kahler potential of the PQ scalar below 
the messenger scale is given by 

c = J d'ez^xM 2 , (12) 

where the wave-function renormalization factor Z& depends on the scale X at the three- 

n 

loop levelo The dominant contribution to the PQ scalar potential comes from rr&(Q = 
k\4>\)\4>\ 2 . In the opposite case, for k\<p\ ^> M mess , the Kahler potential of the X field 
below the PQ scale is given by 

c = J d*ez x (\<t>\)\x\ 2 , (is) 

where the wave-function renormalization factor Zx depends on the scale \(f)\ at the three- 
loop level. Taking them into account, the PQ scalar potential is given by 

Ak 2 o , ,,,, / M„ 



I ^ m softl0| 2l °g I '-rrrr J for *#| < M mess 

IP fl ^fm n (14) 



1 mess 



where m so ft = ((7s/167r 2 )A mess . This drives the PQ scalar away from the origin. In order 
to stabilize the PQ scalar, we introduce the non-renormalizable superpotential 

W=-^r, 15 

M n ~ 2 

where has a PQ charge — n with a cutoff scale M. Then the scalar potential of the PQ 
sector is given by 

0) = v.M) + vu + W + - 2 M 2 ) - ( A Mwk + h - c ) > ( 16 ) 



If there is a mixing between the PQ quarks and messenger fields, the correction arises at the one-loop 
level [35] ■ We do not consider such a case in the following. 
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where V grscv ~ 7n|/ 2 (|0| 2 -|-|0| 2 ) represents the gravity-mediation effect and Am ~ max{m 2 /A r 
m 3/2} comes from the gauge- mediation effect and the gravity-mediation effect, where mo- 
denotes the saxion mass defined later. We choose A M real and positive by the phase 
redefinition of the fields. By minimizing the potential, if the gravity-mediation effect is 
small enough, we find the PQ scale as 



'^2k 2 (M 



fa 



mr 2 sott & V kf a 

—K 2 \F x \ 2 M 2n -Hog 2 ( kfa 



2 n \ ikf mess 



if kf a < M mess 

l/(2n) (17) 

if fc/ a > M mess 



where e = g 4 /(47r 2 ) 3 . Hereafter we take k = 1 and k = 1 for numerical evaluation for 
simplicity, unless otherwise stated. The (ft field also obtains a VEV due to the A-term, as 



n 2fn-2 

Notice that this should be much smaller than f a in order for the calculation so far to 
remain valid. Actually, v/f a ~ Aw/m^ <C 1 is always satisfied. If the X dominantly 
breaks the SUSY, we have a relation 

/^exlO^GeVf^) 173 ^ 173 , (19) 
J VlkeV/ \Mp ' V ; 



for n = 3 and 



for n = 4, if A;/ a 3> M mess . Thus we can have a correct value of the PQ scale for 
M ~ Mp. If the dominant SUSY breaking is carried by another field, we can obtain a 
correct PQ scale for larger gravitino mass. In the large gravitino mass limit, V grav tends 
to dominate the potential and it determines the PQ scale. Fig. @] shows the PQ scale f a 
as a function of M mess for n = 3 (solid) and n = 4 (dashed). The three lines correspond 
to M/Mp = 1, 10 2 , 10 4 from bottom to top for n = 3, and M/M P = 10" 4 , 10" 2 , 1 from 
bottom to top for n = 4. We have taken A mess = 100 TeV and Fx = y/3m 3 / 2 Mp. 

The angular component of the </> around the minimum /„ is regarded as the QCD axion, 
which dynamically solves the strong CP problem^] On the other hand, the fluctuation in 



5 To be more precise, it is a linear combination of the angular components of tf> and <fi that is regarded 
as the axion. In the limit (cf>) ^S> ((f)), however, it mostly consists of the angular component of <j>. 
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the radial direction, which we denote by a, is called the saxion. In the present setup, the 
saxion mass is given by 



p- 1 



It is evaluated as 



for n = 3 and 



= 174 GeV I Io4v) (w)> < 22 > 

fa \* (M P S * 



095 &V lira if!' < 23 ) 



for n — 4. If the X dominantly breaks SUSY, we obtain the saxion mass by using (1191) 
and flUD as 

for n = 3 and 

— v(^) 3/4 (^)" 

for n = 4, if fc/ a ^> M mesa . The axino, a fermonic superpartner of the axion, obtains a 
mass from the operator (TT5|) . Substituting (fl8l) back into ([151) . we obtain the axino mass of 
m„ ~ m a / a/2, taking into account the mixing of 4> and </>Jj The axino also obtains a mass 
radiatively at the three-loop level, but the contribution is suppressed by the three-loop 
factor and the ratio M mess / f a compared with the gaugino mass, which is much smaller 
than the tree-level mass. 

3.2 Cosmology 

Let us discuss the cosmological implications of the present model. In particular, we focus 
on the behavior of saxion. The saxion decays into the axion pair, gauge boson pair, SM 
fermion pair and the Higgs boson pair as long as they are kinematically allowed. The 
decay rates of saxion and are summarized in the Appendix. Since we are interested in 
the saxion mass range m a < 100 GeV, the decay into the SUSY particle pair, including 
the axino pair, is kinematically forbidden. Let us discuss the cosmology in two cases : (i) 
the saxion is initially trapped at the origin, and (ii) the saxion is initially far from the 



6 The 4> and 4> mostly mix with eath other in the mass eigenstate. We call them as the axino and 
denote them by a. 
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Figure 4: The PQ scale as a function of M mess for n = 3 (solid) and n = 4 (dashed). 
The three lines correspond to M/Mp = 1, 10 2 , 10 4 from bottom to top for n = 3, and 
M/Mp = 10 -4 , 10 -2 , 1 from bottom to top for n = 4. 

origin. 

3.2.1 The saxion trapped at the origin 

First let us suppose that the saxion obtains a positive Hubble mass squared in the early 
Universe: V ~ H 2 \tfi\ 2 . It sits at the origin during and after inflation, and the PQ quarks 
are massless there. Therefore, PQ quarks are thermalized and give the thermal mass 
to the PQ scalar <fi. The situation lasts until the temperature drops to the weak scale, 
where the instability of the PQ scalar develops and the PQ symmetry is broken. In such 
a case, the potential energy at the origin dominates the Universe before the PQ phase 
transition, which leads to a short-lasted period of inflation : thermal inflation [551, 13T] . 
(Thermal inflation in the context of PQ symmetry breaking was discussed in Refs. [381 
[39| |4"0~1 |4"H 142] . ) All the dangerous relics such as the gravitino and axino are diluted away. 
The Universe is reheated by the decay of the saxion. For successful reheating, the saxion 
must decay dominantly into fermions or Higgs bosons before BBN begins. Otherwise, 
decay-produced axions contribute too much to the effective number of neutrino species. 
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For example, the saxion may dominantly decay into rf pair for m a ~ 4 GeV, and into 
bb pair for m a ~ 9 GeV. The successful reheating is achieved for such mass ranges. Note 
that there are no axionic domain walls after the phase transition since in this model the 
domain wall number is equal to one. According to the most recent estimate [43j, the PQ 
scale is restricted as f a < 2 x 10 10 GeV taking account of the abundance of the axion 
emitted from the axionic strings and collapsing domain walls if the saxion decays before 
the QCD phase transition. If the saxion decays after the QCD phase transition, the axion 
is diluted by the saxion decay and the constraint is relaxed [44J. 

The field also begins to oscillate around its minimum with an amplitude of v after 
thermal inflation ends. Although it has longer lifetime than the saxion by the factor 
~ (fa/v) 2 (see Appendix), the abundance is suppressed by the factor ~ (v/f a ) 2 . As 
a result, decays before it comes to dominate the Universe. Therefore, it does not 
drastically change the picture as long as <f> decays before BBN. Otherwise, the <f) decay 
can have problematic effects on BBN. 

One should also care about the existence of heavy stable PQ matter, \I/pq, which once 
were in thermal equilibrium during thermal inflation, supporting the PQ scalar at the 
origin. In particular, the mass of the neutral component in \1/pq is severely constrained 
in order for it not to be overabundant [45J. 

Baryon number is also diluted away by thermal inflation. A variant type of the Affleck- 
Dine mechanism may work for generating the baryon asymmetry again after thermal in- 
flation HSl SZl SHI SO]. 

3.2.2 The saxion far from the origin 

Next we consider the case where the saxion is displaced far from the origin because of 
the negative Hubble correction : V ~ — H 2 \(p\ 2 . Then the PQ scalar tracks the minimum 
|0| ~ (i7M n_2 ) 1 ^ n ~ 1 ^ as the Hubble parameter decreases after inflation. In this case the 
saxion relaxes to the true vacuum while the PQ symmetry is never restored. The saxion 
oscillation is induced at H ~ m a , with an amplitude of ~ f a around the minimum. The 
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saxion abundance is then given by 

IT ~ i r " iw^f ~ 2 x irl " GeV (lo^v) (lo4vA (26) 



where e = 1 for T R < y/m a Mp and e ~ y/m^Mp/T-^ for Tr > y/m a Mp. Depending 
on the saxion lifetime and the branching ratio into visible particles, we obtain an upper 
bound on the reheating temperature T R [39]. The saxion decay temperature, T CT , is given 
by 



/ in \ V4 
\7r z g*(T a )J 



'3 G ev(^) 2,3 (^V /6 fern = 3 
VlkeV/ V M / 

* ( ™3/2 \ 7 / 8 /M P \ 5/4 

• ■ <; " w , k :.v) uf' , - " 1 



(27) 

where is the total decay width of the saxion. In this expression we have assumed 
that the saxion dominantly decays into axion pair, and also substituted (1191) and (j2DJ. If 
the saxion decays before BBN, the bound comes from the requirement that the effective 
number of neutrino species, AN e g, must not be much larger than one. However, it does 
not pose a severe constraint. On the other hand, if the saxion decays after BBN, it may 
affect the primordial light element abundances through its hadronic or radiative decay 
processes and hence the saxion abundance is constrained [50J. 

Besides saxion, the coherent oscillation is also induced. As the VEV of decreases, 
the 4> tracks the temporal minimum determined by <p (~ ^4M^ n ~ 2 /|0| n_2 )- At H ~ m a 
where the <fi begins to oscillate around the minimum, the also begins to oscillate around 
the true minimum. The typical oscillation amplitude is estimated to be v given in f ll8p . 
which is much smaller than f a . Although the abundance of is much smaller than the 
saxion coherent oscillation, its lifetime is much longer, as shown in Appendix, and hence 
it is nontrivial whether poses a severe constraint. Actually, as will be shown below, 
can give severer constraint than the saxion depending on their masses. 

The axino (a) is produced by scattering of particles in thermal bath. Since we have 
PQ quarks that couples to the PQ scalar, axinos are produced through the axino-gluon- 
gluino interaction during the reheating [51J. The axino abundance produced by the gluon 
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scattering is proportional to the reheating temperature, and given by 

- 2 x 10-V S (ttJ?—) 2 (77^7) ■ (28) 



In addition, the axino has a tree-level coupling to the Higgs fields and vector-like mat- 
ter. The scattering of Higgs fields and vector-like matter produce axinos as long as the 
temperature is higher than the masses of the Higgs/higgsino and vector- like matter. This 
contribution is roughly given by 



^- 10 " 5 (lo4v)" 2 (^w) 2 ' < 29 > 

where m vec represents the higgsino mass (//) or the mass of vector-like particles, whichever 
is heavier. The axino decays into the gravitino, with the rate 

1 m~ , , , , / m-„ \ 5 / 1 keV \ 2 , . 

r(5->aV>3/ 2 = 5—^—75- ~ (1.1 x 10 sec) {—pttt) • 30 

v ^ i/2! 96vr m\ l2 M 2 p v ' VIGeV/ \m 3/2 J K ' 

Since the axino decays into the gravitino and axion, the bound reads m3/?Ya ^ 4 x 
10- 10 GevEl For relatively small axino mass and large gravitino mass, the axino lifetime 
becomes so long that it dominates the Universe before it decays. Thus, the constraint 
that axions produced by the axino decay do not contribute too much to the iV e g also gives 
upper bound on the reheating temperature. 

Note also that a similar process results in the saxion thermal production and the saxion 
abundance is comparable to the estimates ( 1281) and ( 129]) . Moreover, particles are also 
produced similarly, whose abundance is suppressed by the factor ~ {v/ f a ) 2 . Although 
the abundance of <f) is much smaller than the saxion and axino, its lifetime is much longer 
and hence it is nontrivial whether thermally produced poses a severe constraint. 

We derive cosmological constraints on the present model taking into account all the 
above mentioned contributions : the coherent oscillation of saxion and </>, thermally pro- 
duced saxion, and axino. We follow the methods in Ref. [49J to derive these constraints, 
using decay rates calculated in Appendix. The upper bound on Tr is obtained from the 
requirement that their decay products do not contribute to AN e s and DM abundance too 



7 The gravitino produced by the axino decay has a long free-streaming length and behave as warm/hot 
dark matter rather than the cold dark matter, depending on their masses. In this case the constraint is 
severer. 
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much, do not disturb BBN, do not distort the blackbody spectrum of cosmic microwave 
background, do not yield too much X(7)-ray background. The bound from gravitino 
overproduction [5l] is also considered. See Ref. [49] for details. 

Fig. |5] shows constraints on the reheating temperature Tr as a function of M mess for 
n = 3 and M = 10 4 Mp (top) and n = 4 and M = Mp (bottom). Each line corresponds to 
the bound from the axino, saxion, <p and gravitino. We have also taken Fx = \/?>m^/2Mp. 
A characteristic behavior of the axino bound comes from the fact that, for lower messenger 
scale, the axino becomes light and its lifetime is too long and hence it tends to dominate 
the Universe. For higher messenger scale, the gravitino becomes heavy and its abundance 
coming from the axino decay tends to be too large. Therefore, the bound from axino is 
relatively weak at intermediate messenger scale. The bounds from saxion and <f) shows 
a complicated behavior since their lifetimes significantly change at the threshold for the 
decay into quark/leptons and also the BBN and other constraints significantly depend on 
the lifetime of decaying particles. In the most parameter space, the reheating temperature 
is bounded as Tr < 1 TeV in order to avoid the axino and gravitino overproduction. 
Constraints from the saxion and 4> are less stringent. 

Finally, we mention a constraint on the inflation model. The axion obtains an isocur- 
vature fluctuation during inflation in this case. The magnitude of the CDM isocurvature 
perturbation is given by 



1*1 = = ^M,. = 2 x !0-^ ( *0 " 2 (^] . ,3!) 



a e n c e ttm 1 / 2 vt c e \ io 8 Gev J \ m 



for n = 3 and 



|S ^s?fe- 4x10 tAw^v) [if) ■ (32 > 

for n — 4, where H m f is the Hubble scale during inflation, Q a and fl c denote the density 
parameters of the axion and CDM, respectively, 9 the initial misalignment angle. The 
axion abundance is given by [55j |56j [57] 

Comparing it with the bound \S C \ < 1.4 x 10~ 5 from the WMAP 7 year result combined 
with the baryon acoustic oscillation and the measurement of the Hubble constant [58] , we 
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Figure 5: Constraints on the reheating temperature Tr as a function of M mess for n = 3 
and M = 10 4 Mp (top), and n = 4 and M = M P (bottom). Each line corresponds to the 
bound from the axino, saxion, and gravitino. 
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obtain an upper bound on the inflation scale. If the axion is the dominant component of 
DM, it may be close to the current upper boundB 



4 Conclusions and discussion 

In this paper we have presented a Peccei-Quinn extended GMSB model in which the extra 
vector-like matters are introduced in order to explain the 125 GeV Higgs boson. The PQ 
symmetry controls the size of /x-term as well as the masses of the vector-like matters. 
These parameters have correct size for the PQ symmetry breaking scale of 10 10 -10 12 GeV. 
Thus the strong CP problem as well as the SUSY CP problem are solved in our model. 
Fortunately, the PQ symmetry forbids the unwanted coupling between the vector-like 
matter and the down type Higgs, which would otherwise make negative contribution to 
the Higgs boson mass and induce a CP violating phase. There are parameter regions 
where the muon g — 2 anomaly is explained by the SUSY contribution consistently with 
the 125 GeV Higgs boson. We have also constructed a method to stabilize the PQ scalar 
and shown that the correct PQ scale can be obtained for natural parameter choices. We 
have discussed cosmological effects of the saxion and axino, and derived upper bound on 
the reheating temperature. 

Some notes are in order. Although we have focused on the GMSB model, a similar 
model can also be applied to the gravity-mediation models such as the CMSSM. The Higgs 
mass as well as the muon g — 2 can be explained simultaneously. Although there are no 
messenger fields, the stabilization of the PQ scalar is achieved by the balance between the 
negative gravity-mediated mass term V ~ — m^ 2 \<f)\ 2 and the non-renormalizable potential 
arising from Eq. ( TT5|) . As another method, we may simply introduce the superpotential 
as 

W = K'S(^-f 2 a ), (34) 

where S is a singlet field. In this case, the positive gravity-mediated mass term V ~ 
m 3/2(l^l 2 + l^l 2 ) stabilizes the PQ scalars at \<fi\ ~ \<f)\ ~ f a . This may also be consistent 
with the hybrid inflation model of Ref. |60j . 

8 If the misalignment angle 9 is small while keeping the magnitude of S c , there can be large non- 
Gaussianity in the CDM isocurvature perturbation |59j . 
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A Decay of saxion and (j) 

In this appendix we list up the decay modes of the saxion and <fi. First, we consider the 
mixing between PQ scalars, (ft and (p. Let us expand them as 

= fa + -j=(a + ia), (35) 

4> = v H — -=(a + ia). (36) 
V2 

They couple in the superpotenial as W = n 0/M™ -2 . By noting the relation (fl8l) . the 
scalar potential ( IT6l) is expanded as 



V = + Wa + -^Ua 2 + -^Laa 2 + ( " = ( l<xa 2 + iaa 2 ) ,(37) 



2 CT 2 2a/2/ 2^ v^n V/a 

where a = (a, a), a = (a, a) and 

^ = fn- 2 TL ^1^): (38) 



KA i fnA M m & (v/ f a ) A M m- a . , 

a 1 A M m f (l/n)A M m s (/ a /^) ' ' 



where 



2 1 d 2 ^(|0|) 

777 = — 

2 ;)LJL> 

j.2(n-l) 
2 2 J" 

m= = n 



M 2 ("- 2 ) 



(40) 



]Vf2(n-2) • 

As for the CP-odd parts, there is a massless mode since det(A^ 2 ) = 0, as is expected from 
the PQ symmetry. In what follows, we approximate as m a ~ v^n/;- 1 /!"- 2 . Then we 
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obtain 




(42) 



(41) 



For f a ^> v, the mass eigenstates are given by 



nv _ 

a ~ a — —a, 



(43) 



J a 

nv 



(44) 



a ~ a + -7- a, 



where a is the massless Goldstone mode, which is regarded as the axion, while a has a 



where the saxion a, which mostly consists of a, has a mass of ttv and a has a mass of 
m a /V2. 

A.l Saxion decay 

Here we summarize the saxion decay mode. Since the saxion is at most around (9(100) GeV 
in our model, we neglect the decay into SUSY particles. We also ignore the mixing of a 
and a in this subsection, since it does not affect the saxion decay rate summarized in the 
following as long as the mixing angle is small. 

A. 1.1 Decay into axions 

From the Lagrangian (IBTj) . the saxion decay rate into axions is calculated as_ 



mass of m a / a/2. In the most part of this paper, we have not distinguished a and a since 
the mixing angle is small. 

As for the CP-even part, the mass eigenstates read 




(45) 



(46) 




(47) 



9 If the PQ scalar is expanded as <fi = /aexp [(a + ia)/^/2j 
kinetic term C = \d(j)\ 2 . 



we find the same decay rate from the 
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Numerically, it is evaluated as 

r(a^2a)^(lxlO- 2 sec)- 1 ( T ^j ^ - j . (48) 

A. 1.2 Decay into gluons 

if the saxion is heavier than ~ 1 GeV, it can decay into the gluon pair, which hadronize 
and results in the production of energetic particles. Integrating out the PQ quarks yields 
the following couplings 

This induces the decay into gluons as 

32vr3 ft 



V{a^2g)^^^. (50) 



A. 1.3 Decay into fermions 

Let us expand the Higgs and PQ scalar as 



H u = v u + ^=(h u + ia u ), (51) 
H° d = v d + -f=(h d + ia d ). (52) 



They have a coupling in the superpotential as 



w = W~ iHuIId - (53) 

We have focused on the case of £ = 2 in this paper, but here we do not fix it. The /i-term 
is generated through the VEV of as /x = A/^/Mp" 1 . Let us focus on the CP-even parts. 
First, by diagonalizing the h u and h d , we obtain the light and heavy Higgs bosons as 

h = h u cos a — h d sin a, (54) 
H = h u sin a + h d cos a. (55) 
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They mix with the saxion a in the mass eigenstates. Writing the mass eigenstates as h 
and H, we find 



h ~ h 
H ~ H + 



ifXV 



fa(m 2 h -m 2 a ) 



[B cos(a + 0) + 2/1 sin(o; — 0)] cr, 
[-B sin(a + 0) + 2/j, cos(a — 0)\ a, 



(56) 
(57) 



f a (m 2 H - ml) 

where v = (v 2 + v^) 1 ^ 2 . From these mixings, we can calculate the saxion decay rate into 
the SM fermion pair. For up(down)-type quarks, we obtain 



r(<7 ->■ f u (d)fu(d)) 



3£ 2 m a m 2 f c u{d) n 



.2 \ 2 



16tt fl \m 2 h - ml 
where m/ denotes the final state fermion mass and 



1 - 



4m 



2x 3/2 

/ 



mi 



(58) 



cos a 



c d 



■"U, ■ n 

sin p 
sin a 



cos (3 



— cos(o; + 0) + 2 sin(a — /3) 



+ — 



2 2 



cos(a + /3) + 2 sin(a — /?) 



+ 



sin p mjj — m* 
cos a ml — ml 
cos j3 m 2 H — ml 



— sin(a + 0) — 2 cos(a — 0) 



B 
// 



sin(o; + 0) — 2 cos(a — /3) 



For the decay into charged leptons, we obtain 



r(<r ->• ZZ) 



£ 2 m a m 2 f c d ji 



1 



4m 



2\ 3/2 



mi 



(61) 



167T /2 ^_ m 2^ 

where m/ is the final state lepton mass. 

For later convenience, we describe the mixing of CP-odd parts. They form mass 
eigenstates as 



G = a d cos j3 — a u sin/3, 

Q iv sin(2/3) 
a H = a d sm p + a u cos p H — a, 

fosin(2/3). . fl 
a = a — (a d smp + a u cosp), 



(62) 
(63) 

(64) 



where G corresponds to the Goldstone boson eaten by the Z-boson, au is identified as 
the CP-odd Higgs boson, and a is the massless mode in association with the spontaneous 
PQ symmetry breaking. 
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A. 1.4 Decay into Higgs bosons 

If the saxion is heavier than the twice the (lightest) Higgs boson mass, the saxion can 
decay into the a pair of the Higgs boson. The decay rate is given by 



n 2 ml ( a \ 4 / B sin(2a)\ 2 / 4m 2 h \ 1/2 
A. 2 (f) decay 

Next let us estimate the decay rate of 0. It decays into the axion pair, gauge boson pair 
and the SM fermions through the 0-0 mixing. Roughly speaking, the decay rates of a 
and a are suppressed by the mixing factor ~ {nv/f a ) 2 compared with the saxion decay 
rate, since does not have direct couplings to light particles. 

A. 2.1 Decay into axions 

The a decay rate into two axions can be read from fl37j) . taking into account the mixing 
between a and a, as 

r ( ^2 a)= <» 2 + 2 "- 2 > 2 ^y. (66) 



6*r fl \U 
A. 2. 2 Decay into gluons 

From Eq. f l4"9"|) . the a(a) decay rate into the gluon pair is obtained as 
and 

(S ■ (68) 

A. 2. 3 Decay into fermions 

The a decays into fermions through the mixing of a— a. From Eq. (jlp]l . we find the decay 
rate into a fermion pair as 

r(^//^p^)V^//), (69) 
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after the m a in the formula is replaced with m„. Similarly, we can find the decay rate of 
a into an up(down)-type quark pair from the mixing given in (jUJ), as 

n^ fM = -^i (l-/) , (70) 



where 



, (tan/3) 1 , tan/3 

° u = tan/3 + (tan/3)" 1 ' ° d = tan/3 + (tan/3)- 1 ' ( ^ 



The decay rate into a charged lepton pair is given by 



£ 2 m a m 2 ( c' d nv\ 2 ( Amf^ 



na^U) = —^k L [^) (1--^) • (72) 



B Neutrino mass 



All gauge-invariant and R-parity conserving operators up to dimension four in the super- 
potential are listed in Table. [20 It is seen that the neutrino mass operator LH U LH U is 
allowed for qu = —4/5. Actually, the operator 

W = ^M N NN + y N NLH u , (73) 

is allowed for qu = —4/5 if the right-handed neutrino has a zero PQ charge. In this case, 
the higher dimensional operators 

W = W ( QQQ ' Hd + QU'EH d + QUE'H d ) (74) 

are also allowed and the vector-like particles can decay (M* is the cut off scale). The 
Kahler potential like 

K= W ( QU ' L] + UE ' D1 + E ' ]HdHd + E ' H d H D + h - C - ( 75 ) 

are allowed, which also induce the decay of vector-like matter. Note that the life time of 
the proton is sufficiently long, since its decay width is suppressed by M*~ 4 . 
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